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We present the zero-temperature phase diagram of the one-dimensional t2g-orbital Hubbard 
model, obtained using the density-matrix renormalization group and Lanczos techniques. Emphasis 
is given to the case for the electron density n=5 corresponding to five electrons per site, of relevance 
for some Co-based compounds. However, several other cases for electron densities between n=3 and 
6 are also studied. At n—5, our results indicate a first-order transition between a paramagnetic (PM) 
insulator phase and a fully-polarized ferromagnetic (FM) state by tuning the Hund's coupling. The 
results also suggest a transition from the n=5 PM insulator phase to a metallic regime by changing 
the electron density, either via hole or electron doping. The behavior of the spin, charge, and orbital 
correlation functions in the FM and PM states are also described in the text and discussed. The 
robustness of these two states varying parameters suggests that they may be of relevance in more 
realistic higher dimensional systems as well. 

PACS numbers: 75.30.kz, 71.10.Fd, 75.50.Cc, 75.10.-b 



I. INTRODUCTION 

The study of the exotic properties of cobalt oxides is an 
area of investigations that is currently attracting consid- 
erable attention in the research field of condensed mat- 
ter physics. Among the main reasons for this wide ef- 
fort, the recent discovery of superconductivity in layered 
two-dimensional triangular lattices of Co atoms with the 
composition NaxCo02 has certainly triggered a rapid in- 
crease of research activities on cobalt oxides. This mate- 
rial becomes superconducting after H2O is intercalated^ 
opening an exciting area of investigations. 

The experimentally unveiled phase diagram of this 
compound varying the Na composition has revealed the 
existence of several other competing tendencies: Charge 
ordered as well as magnetic states are stabilized, in ad- 
dition to superconductivity^ In the related compound 
(Ca2Co03)(Co02), a spin incommensurate spin-density- 
wave has been recently reported.'^ The existence of such 
a rich phase diagram is a characteristic of strongly corre- 
lated electron systems, where complex behavior typically 
emerges due to the presence of competing states that 
have similar energies but vastly different transport and 
magnetic properties 1^ 

Additional motivation for the study of Co-oxides arises 
from recent experimental studies of hole-doped cobaltites 
in the perovskite form such as Lai-xSrxCoOa, where 
clear tendencies toward phase separation between ferro- 
magnetic (FM) metallic and paramagnetic (PM) insulat- 
ing regions have been found. ^ This establishes an intrigu- 
ing qualitative connection between Co-oxides and the fa- 
mous manganites that exhibit the colossal magnetore- 
sistance, effect widely believed to originate in an analo- 
gous mixed-phase tendency exhibited by Mn-oxidesA In 
fact, a large magnetoresistance in some cobaltites has 



also been observed, and its origin appears related with 
phase competition.^ As a consequence, establishing the 
dominant ground-state tendencies of simple models for 
cobaltites is important to envision the possible phase 
mixtures that may lead to exotic behavior. 

As a third motivation for the study of cobaltites, it is 
known that some Co-based compounds have interesting 
thermoelectric properties. In particular, a huge thermo- 
electric power has been recently discovered in NaCo204 
by Terasaki et al.^ opening another area of investiga- 
tions, with a focus on thermoelectric materials mainly 
with the purpose of industrial applications. 

For all these reasons, the theoretical study of mod- 
els for Co-oxides is timely and needed in order to guide 
further experimental developments. Ab-initio calcula- 
tions have already provided important information in 
this context)^ and the inclusion of many-body effects is 
the natural next step. Previous theoretical studies of 
Co-based systems including Coulombic repulsion have 
mainly focused on triangular lattices. In this context, 
recent Monte Carlo investigations unveiled the presence 
of magnetic correlations. Fluctuation-exchange approx- 
imations also revealed tendencies toward ferromagnetism 
and possible triplet-pairing instabilities in a multiorbital 
modeliii Several approximate studies of t-Ji^ and single- 
band Hubbard models'-' have also been presented. 

To understand the behavior of complex oxides, it is 
of particular importance the analysis of the many possi- 
ble tendencies in the ground state, namely the study of 
the various competing states stabilized as electron den- 
sity and coupling are modified. Unfortunately, this task 
is difficult due to a lack of reliable unbiased analytical 
techniques. For this reason, the first effort toward a de- 
tailed numerical analysis of models for cobaltites is pre- 
sented here. Instead of directly emphasizing the triangu- 
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lar lattice with approximate techniques or exactly study- 
ing small systems, we have preferred to perform a sys- 
tematic study of a one-dimensional multiorbital Hamilto- 
nian, exploring in detail the coupling and electron density 
parameter space, and using computationally exact tech- 
niques. This level of accuracy is achieved through the use 
of reliable methods such as the density-matrix renormal- 
ization group (DMRG)^"' and the Lanczos techniqueii^ 
We envision this effort as a first step toward a system- 
atic computational analysis of more complicated quasi- 
two-dimensional triangular-lattice systems. 

The paper is organized as follows. In Sec. II, the 
multiorbital model is introduced and many-body com- 
putational techniques used here are briefly discussed. In 
Sec. Ill, the main results are presented. These results 
are organized based on the observable studied: First, the 
n=5 phase diagram is discussed, where n denotes the 
number of electron per site. Then, the spin correlations 
are presented at several values of n. This is followed by 
the charge and orbital correlations. Finally, conclusions 
are presented in Sec. IV. The main result of the paper is 
the clear dominance of two rather different ground states: 
(1) a FM state and (2) a PM state with short-range cor- 
relations. Both are very robust varying couplings and 
densities. Their higher-dimensional versions may be of 
relevance for present and future Co-oxide experiments. 

II. MODEL AND TECHNIQUE 

In the investigation reported in this manuscript, we 
consider a three-orbital Hubbard model, defined on a one- 
dimensional chain along the x-axis with L sites. The 
three orbitals represent the t2g orbitals of relevance for 
cobaltites. The model is given by 

^ = - E ^7,7' (47.'^,+i,y.+H.c.) 

E ^],70-'^l,7'(T''^j,7<^''^j,7'o- 

+ y E 4,7't4,70-''^J^7''^"^J,7'<^' (1) 
J,o-#o-',75^7' 

where the index j denotes the site of the chain, 7 indi- 
cates the orbitals xy, yz^ and zx, and a is the spin pro- 
jection along the 2-axis. The rest of the notation is stan- 
dard. The hopping amplitudes are ixi/.a;j/=^za:,zK=^=l, 
and zero for the other cases. These simple values for the 
hopping amplitudes can be easily derived from the over- 
lap ofdxyj dyz, and dzx orbitals between nearest- neighbor 
sites along the x-axis. The interaction parameters U, U' , 
J, and J' are the standard ones for multiorbital Hamilto- 
nians, and a detailed description can be found in Ref. G. 
These couplings are not independent, but they satisfy the 
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FIG. 1: Ground-state phase diagram for the one-dimensional 
three-orbital Hubbard model, using a 6-site chain and work- 
ing at electron density n=5. FM and PMI denote the regions 
with ferromagnetism and paramagnetism (insulator), respec- 
tively. We also present a schematic picture of the electron 
configurations. AFO indicates the staggered population of 
orbitals in the FM state. The reader should consult the text 
for more details, as well as the next figure. 

well-known relations J'=J and U=U'+2J, due to the re- 
ality of the wave function and the rotational symmetry 
in the orbital space. 

We investigate the model described above mainly using 
the DMRG technique with open boundary conditionsii^ 
The finite-size algorithm is employed for sizes up to 
L=48, keeping up to m—350 states per block. The trun- 
cation errors are kept around 10~^ or smaller. The center 
blocks in our DMRG procedure are composed of 64 states 
due to the three orbitals. Note that, for instance, the t- 
J model has only 3 states in these center blocks. As a 
consequence, keeping to=350 states per block in the t2g- 
orbital Hubbard model is analogous to keeping m^7000 
states per block in the t-J model. 

Although in related investigations specific values of U 
and J for the triangular Co-oxides were discussed, here 
we prefer to vary independently these couplings analyz- 
ing the possible ground states that are stabilized by this 
procedure. In fact, the ratio U/J may change among 
the many interesting Co-oxides, and in addition, it is im- 
portant to classify the states that could be stabilized by 
the proper isovalent chemical doping, external fields, or 
perturbations. 



III. RESULTS 
A. Phase diagram for density n=5 

In Fig. 1, the ground-state phase diagram J versus 
Uctf—U'—J is presented. For large J, a fully polarized 
FM phase is obtained, while for small J, a PM regime 
is found. This PM phase is insulating, as shown below. 
Note that the phase diagram is obtained by comparing 
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FIG. 2: States with the largest weight in the ground state of 
a 4-site chain solved exactly. Note that each state has 8-fold 
degeneracy. At J—0, these three states have the same weight. 
On the other hand, for nonzero J, the state (a) (and its 8- 
fold degenerate states) has the largest weight, with a spin 
(orbital) structure factor peaked at 7r/2 (tt). The states (b) 
and (c) (each one also with degeneracy 8) have the second- 
and third-largest weights, respectively, for nonzero J. 



the energies for different sectors of the z-projection of 
the total spin, S^^^^i, mainly using a system of size L—6. 
Other values of L are also studied, and it is observed that 
for L—A, 6, 8, and 10, in the PM regime the ground state 
has total spin 0, 1,0, and 1, respectively, for a large set of 
couplings investigated. As a consequence, it is reasonable 
to assume that the transition line separates states with 
the minimum and maximum total spin, without interme- 
diate partially polarized regimes. The phase diagram we 
have found has similarities with that already reported by 
two of the authors at density n=4,— in the context of 
spin-1 chains. 

As described later, our results for the spin-spin, charge- 
charge, and orbital correlations suggest, roughly, an elec- 
tron distribution at short distances, as schematically pre- 
sented in Fig. 1. The electron configuration in the FM 
phase is quite simple: 5 electrons per site, with a po- 
larized net spin 1/2 and antiferro-orbital (AFO) correla- 
tions. The existence of FM correlations is a direct con- 
sequence of the multiorbital nature of the model and the 
robust value of J in the FM regime. 

On the other hand, a more complex electron config- 
uration emerges in the PM phase. The meaning of the 
full circles in the inset of Fig. 1 for the PM phase is to 
denote either a spin up or a spin down. Note, however, 
that quantum fluctuations are strong and the configura- 
tion shown in Fig. 1 is just a guidance. To obtain insight 
into the ground-state wave function, it is useful to con- 
sider the case of a four-site chain, where results can be 



obtained exactly by using the Lanczos method. In the 
strong-coupling limit Ucs ^ J ^ 1 (or, more precisely, 
1/{U' ~ J) ^ 1), it is found that the most important 
portion of the ground-state wave function is expressed as 
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where the sum is taken over the permutation of the four 
spinors and np is the number of permutation we have to 
perform to recover the original configuration. Namely, 
the electron configuration presented in the PM phase 
of Fig. 1 should be regarded as the equivalent of the 4 
spinors contained in lip). Note that this is not a rigid con- 
figuration, but all permutations are equally important at 
small J. In particular, all the 24 states have the same 
weight in the ground state at J=0, while at finite J, the 
24 states are split into three classes with 8 states for each, 
as shown in Fig. 2. Note that each of these classes lead to 
a distinct peak position in the spin and orbital structure 
factors. When the peak positions in these channels are 
denoted by gspin and gorbitai, the class (a) has gspi,i=7r/2 
and gorbitai=7r, class (b) gspin=7r/2 and qorbitai=7r/2, and 
class (c) <?spin=7r and (j'orbitai=7r/2. The yz orbital is fully 
occupied due to the one dimensionality of the system that 
prevents the movement of electrons in this orbital due to 
a vanishing hopping. For the active xy and zx orbitals, 
the electrons are not distributed in a rigid charge-ordered 
pattern, but instead the density is to an excellent approx- 
imation equal to 1.5 at every site. Note that a similar 
representation of the ground-state wave function for four 
sites has been found for the SU (4) spin-orbital model.— 
As discussed in more detail below, these two models are 
related to each other. 



B. Spin correlations at several densities 

To understand more quantitatively the magnetic order 
present in the PM phase, it is useful to measure the spin- 
spin correlation function, defined as 



1 



Cspin{l) — -TJ i^i^j) ' 



M ^ 



(3) 



where Sl=^^{pi^^ — pi^i)/2 is the z-projection of the 
total spin at each site and M is the number of site pairs 
(i,j) satisfying /=|i—j|. We average over all pairs of sites 
separated by distance in order to minimize boundary 
effects. In Fig. 3(a), Cspin(0 is shown for the PM phase. 
It is found that the numerical data of Cspin(0 ^re well 
reproduced by the function 
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as shown by the dashed curve. Note that a and h are 
appropriate fitting parameters. The result indicates that 
the spin-spin correlation function has a four-site period- 
icity and decays as a power law with critical exponent 
3/2. In the inset of Fig. 3(a), we also present the Fourier 
transform of the spin-spin correlation function, 

S{q)^\Y.^"'^'''^{S',Sl), (5) 

for L=16 and L=48. As observed in this figure, finite- 
size effects appear to be very small. Here we clearly find 
a peak in S{q) at g=7r/2, corresponding to the four-site 
periodicity of the spin-spin correlation function. 

Let us here discuss the physical meaning of the four- 
site periodicity. Since the yz orbital is fully occupied 
in our studies, the t2g-orbital Hubbard model can be re- 
garded as a two-orbital Hubbard model composed only of 
xy and zx orbitals. Note that for this two orbital model, 
the hopping amplitudes are symmetric and there is no 
off-diagonal elements. Moreover, when J=0, in the two 
orbital model, there exists an extra SU{A) symmetry in- 
volving both spin and orbital degrees of freedom. In such 
a case, the effective Hamiltonian in the strong-coupling 
limit is given by the SU{A) spin-orbital model, which 
has been investigated intensively in recent years >i2ii2i2£ 
Concerning a less symmetric case than S't/(4), the effect 
of J has also been discussed fSLSS where anisotropic ex- 
change interactions arise in the orbital part. Note that 
the spin-spin correlation function, presented in Fig. 3(a), 
is found to have a four-site periodicity and decay as a 
power law with critical exponent 3/2 These results are 
consistent with previous analytical work^ and numerical 
analysiaiSiSi for the SU{4) spin-orbital model. 

Here we stress that the spin-spin correlation functions 
for n=5 clearly present distinct behavior from the results 
already reported at n=4fi& where an exponential decay 
has been observed, as depicted as a linear-log plot in 
Fig. 3(b). The result indicates a gapless spin-excitation 
spectrum for n—5, with power-law decaying correlations, 
in contrast to a gapfuU behavior for n—4. Note that, for 
better comparison, we have normalized Cgpin in such a 
way that the correlations are the same at distance one. 
We have eliminated the odd sites for n=5, since the re- 
sults there are close to zero (see Fig. 3(a)). Note also that 
working with m=350, it is difficult to reach good accu- 
racy for the correlations at large distances, since they are 
very small. For this reason, in our results we present only 
the first 19 sites 

In Fig. 3(b), we also show the spin-spin correlation 
function for n=3, reported here for the first time to our 
knowledge. In the case of n— 3, it is naively expected that 
the local spin S—3/2 is formed at each site. By analogy 
with the half-odd-integer-spin antiferromagnetic Heisen- 
berg chains, we expect the power-law decay of the spin- 
spin correlation function and a gapless spin-excitation 
spectrum as well^ On the contrary, we can observe in 
Fig. 3(b) that the spin-spin correlation function shows an 
exponential decay similar to the case for the integer-spin 
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FIG. 3; (a) The spin-spin correlation function Cspin(j) vs. 
j for L—48 and density n=5. The dashed line indicates a fit 
using Eq. |0J . The inset shows the spin structure factor S (q) 
for L—16 and L—48. (b) The linear-log plot of the module of 
the spin-spin correlation |Cspin(OI corresponding to densities 
n=3, 4, and 5 with L=48, as well as the fit used in (a). For 
details, see the main text, (c) Spin structure factor S (q) for 
several densities n, and using L=16. The arrows indicate the 
peak positions. In all plots i7cft=10 and J—1, as indicated. 
Inset shows S{q) for n=4 and 5. 



chains. To understand this peculiar behavior, it is nec- 
essary to take into account the effect of t2g orbitals. As 
mentioned above, electrons in the yz orbital cannot hop, 
while electrons in the xy and zx orbitals move to adjacent 
sites with the same amplitude. Then, it is expected that 
only electrons in the xy and zx orbitals contribute to the 
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exchange interaction, and the n=3 system is regarded as 
an effective 5=1 chain, leading to the exponential decay- 
ing spin-spin correlation function. 

In Fig. 3(c), S{q) is shown for several densities. Since 
the finite-size effects seem to be small, we consider L=16. 
As observed in these studies, the results suggest that 
the peak position changes linearly with the electron den- 
sity as q={6 — n)7r/2 (mod tt). Note that this peak is 
clearly robust for n=4, as shown in the inset of Fig. 3(c), 
and substantially decreases its intensity by increasing the 
density n. 

It is important to remark that the inset of Fig. 3(a) 
is very similar to the results found by Ogata and Shiba 
in their study of the one-dimensional Hubbard model at 
quarter-filling and U=oo (see Fig. 9 of Ref . l27f) . Clearly, 
in the model studied in this paper, the electrons in the 
two bands with a nonzero hopping behave like one-band 
models with a strong on-site repulsion, at least from 
the perspective of the spin correlations. Note, however, 
that these two one-band models are connected via the 
Coulombic repulsion which, as discussed below, will open 
a gap in the spectrum of charge excitations. 



C. Charge correlations at several densities 

To investigate the charge excitations, it is useful to 
measure the charge gap, defined as A=E{Ne -I- 2) -|- 
E{Ne - 2) - 2E{Ne), where E{Ne) denotes the lowest en- 
ergy in the subspace with the total number of electrons 
N^. In Fig. 4(a), the charge gap is shown as a function 
of 1/L at densities n—A and 5, for particular values of 
C/eft and J. Clearly, at these densities the charge gap 
extrapolates to a nonzero value in the thermodynamic 
limit, indicating that the system is an insulator. On the 
other hand, as shown in Fig. 4(b), for non-integer elec- 
tron densities, the charge gaps seem to extrapolate to 
zero in the thermodynamic limit, suggesting a metallic 
behavior. These results indicate that a transition from 
an insulating phase to a metallic regime is obtained by 
changing the density away from n=5. 

In Fig. 4(c), the charge gap for the density 7i=5 and 
L=12 is presented. It appears that U' is the main driver 
of the system into an insulating phase. On the order 
hand, the Hund's coupling J has the opposite effect: As 
observed in Fig. 4(d), by increasing J the charge gap 
decreases. Note that U' plays a role similar to that of 
the nearest-neighbor Coulomb repulsion V in the two- 
leg ladder extended Hubbard model (with the two legs 
playing the role of the two orbitals in our model) . In the 
ladder case, it is known that V drives the system to an 
insulator at quarter-filling 

We have also investigated the charge structure factor, 
defined as 

^'^'^'(9) = ^E^''^'''^ (W(j,fc) +7V^''^(j, fc)) , 

(6) 
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FIG. 4: (a) The charge gap A vs. l/L at particular values 
of UcH and J, and densities n=4 and n—5. (b) Same as (a) 
but for non-integer densities, and Ucft=10 and J—1. (c) and 

(d) denote the charge gap for density n=5 and L=12. (c) 
contains A vs. U' at J=l, while (d) shows A as a function of 
J at U'^11. 

where N"''^ {j,k)={5n^{j)5n^i(k)) and 5n^{j)=n^{j) — 
(%(j))- In a periodic system TV'"^ {j,k)=N'^ '''{j,k). 
However, with open boundary conditions, as used in our 
investigation, this is not valid any more, due to the pres- 
ence of Friedel oscillations. Using the definition discussed 
above, N'^''^ {q — 0) is always zero. In our calculations, 
we obtained TVT'^ (g = 0)<10^^, indicating that we have 
retained enough states in the truncation process to sat- 
isfy this constraint. 

The best indication of a true long-range-order (LRO) 
can be obtained by the system-size dependence of 
N''-'^'{q). liN'i''<' [q*)/L constant as L ^ oo, at some 
particular q* ^ a true LRO characterized by q* is present. 
Carrying out this analysis, we have found no evidence of 
LRO in the charge sector of n=5. In Figs. 5(a), (c), and 

(e) , typical examples of the charge structure factor for 
the FM and PM phases at density n=5 are presented. In 
the FM phase, we are able to explore very large system 
sizes, since we can measure the correlations in the sector 
of 5'jQj3^j=max, with a much smaller Hilbert space than 
for the PM phase. Although we did not find LRO, the 
behavior of the structure factor suggests that in the FM 
phase the charge-charge correlation presents a quasi-LRO 
due to the presence of a robust peak at q='K. In fact, in 
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FIG. 5: The charge structure factor A^"^'"^ (q) and the charge- 
charge correlation function C{j), at density n=5. (a) and (b) 
are for [/off =10, J=10 and I/=64. This corresponds to the FM 
regime of Fig. 1. The dashed line is a fit using the function 
acos{'Kj)/j with an appropriate fitting parameter a. (c) and 
(d) are for UcS=W, J—1, and L—32. This is in the PM regime 
of Fig.l. (e) and (f) are the same as (c) and (d), respectively, 
but for J=5. 



the charge-charge correlation function, defined as 

^"(0 = ^ E iSn.A^)Sn.Aj)) , (7) 

we observe a slow power-law decay, as shown in Fig. 5(b). 
This correlation osciUates as cos(7rj)/j, as indicated by 
the dotted curve in Fig. 5(b). The DMRG data agree 
very nicely with a fit using this function. These strong 



0.3 



Q-Hn=4.75 

n=5 
^n=5.25 
<«>n=5.5 




FIG. 6: The charge structure factor N'"'^''''^ (q) for several 
densities and using UcB~10, J=1, and L—16. The arrows 
indicate the cusp positions. 



charge oscillations suggest that the system may develop 
LRO rapidly, when a coupling in the direction perpendic- 
ular to the chains is introduced. Then, spin-FM charge- 
ordered states should be seriously considered as a pos- 
sibility for Co-oxide materials, although more detailed 
calculations are needed to confirm this speculation. 

Note also that the negative values of iV^^'^^(g = it) 
suggest an alternation of charge occupation between the 
zx and xy orbitals, as in the schematic representation in 
Fig. 1 (FM phase). Indeed, as discussed later in more 
detail, there is quasi-long-range AFO order. A similar 
result has already been found in the FM phase for the 
density n=Ai^ 

On the other hand, in the PM phase, N'^'f [q) does 
not present a peak as sharp as for the FM phase, as 
shown in Fig. 5(c). In fact, the magnitude of the charge 
correlations is drastically different between the PM and 
FM phases, as can be seen from the absolute values of 
these correlations in the vertical axes of Figs. 5(b) and 
(d). Also note that the appearance of the cusp at g=7r/2 
is related to the four-site periodicity of the correlation 
C{1), as shown in Fig. 5(d). 

Our results also suggest that the charges behave dif- 
ferently in two distinct regimes in the PM phase. At 
small J, the correlation C{1) presents a four-sites peri- 
odicity, while for larger J, only a two-site periodicity is 
found, as observed in Fig. 5(f). In addition, the cusp 
of N^^'^^(q) present in the small- J regime disappears at 
larger J (Fig. 5(e)), apparently continuously. We have 
also observed that at small J, the position of the cusp 
changes with the electron density in a similar way as 
S{q), as shown in Fig. 6. 



D. Orbital correlations at n=5 

Consider now the possibility of orbital order. In the 
PM phase and for n=5, we have found that the xy and 
zx orbitals are those of relevance, since the yz orbitals 
are fully occupied. Note that in the PM phase and at 
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FIG. 7; (a) The orbital structure factor T{q) versus momen- 
tum for (7off=10, J=10, and L=64 with 9i=e^0. (b) The 
orbital correlation Corbitai(j) for the same parameters as used 
in (a). The dashed line is a fit using the function acos{nj)/j. 
(c) and (d) are the same as (a) and (b), but for UcS=W, J=l, 
and L—32. (e) and (f) contain the correlations CspinO) and 
CorbitaiO) for UcH~W, J—O, and L—32. All the results are 
for the density n=5. 



n—A, the orbital degree of freedom becomes inactive due 
to the ferro-orbital order ii& Then, here we take the pseu- 
dospin representation for the xy and zx orbitals, and 
measure the orbital correlations to determine the orbital 
structure. For this purpose, we introduce an angle dj 
to characterize the orbital shape at each site. Using the 
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FIG. 8: (a) The size dependence of the orbital structure 
factor T{q) at q—n with 6i—6—0, at density n=5. (b) T{q) 
vs. 9 for particular values of q. 



angle 9j , we define the phase-dressed operators as 

The optimal orbitals, a and b, are determined so as to 
maximize the orbital structure factor, defined by 



ilU-k) /rpz 



(9) 



where T^U)=EAfla.aU.a..^ - fib,Mb,<r)l^- 

Let us first focus on the case 9i=9=0. In Figs. 7(a) 
and (c), typical examples of the orbital structure factor 
in the FM and PM phases at density n=5 are presented. 
Note that these results are similar to those of the charge 
structure factor shown in Figs. 5(a) and (c), as previ- 
ously anticipated. Also, as shown in Figs. 7(b) and (d), 
concerning the orbital correlation function defined as 



Co 



■bital 



(0 



1 

M 



J2 {T^^)T^J)), 

H-3\=i 



(10) 



we find the same form as C{1), as observed in Figs. 5(b) 
and (d). In the FM phase, as shown in Fig. 7(b), 
Corbitai(0 dccays as cos(7rj)/j, which is the signature of 
quasi-long-range AFO. On the other hand, in the PM 
phase, we observe a four-site periodicity of Corbitai(j) as 
well as that of Cspin(i), while the peak position of T{q) 
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is at q='K for /7off=10 and J=l. Note that the spin-spin 
correlation function shows the four-site periodicity and 
S{q) has the peak at q—TT/2 for Ucf{—1Q and J=l, as 
shown in Fig. 3(a). 

To clarify the similarity between the two-orbital model 
composed of the xy and zx orbitals and the SU{A) spin- 
orbital model, we investigate Cspin(0 ^-nd Coibitai(0 fo'" 
the present t2g model, at [/cff =10 and J=0. As shown in 
Figs. 7(e) and (f), it is clearly observed that Corbitai(0 
and Cspin(0 present exactly the same behavior with a 
four-site periodicity, due to the presence of the SU{4) 
symmetry at J—Q. When we include the effect of J, 
the spin and orbital degrees of freedom are no longer 
equivalent, but we can observe the four-site periodicity 
in both Coibitai(0 ^-i^d Cspin(0 due to the influence of the 
SU{A) symmetry at J—Q. Thus, the short-range orbital 
correlation for small J originates in the SU{A) singlet at 
J=0. 

It should be mentioned that there is no indication of 
orbital LRO in the PM phase, since T(7r) converges to 
a finite value in the thermodynamic limit, as shown in 
Fig. 8(a). On the other hand, although we have found no 
signature of orbital order between the xy and zx orbitals 
through the orbital structure factor T{q) for 6i=6=0, a 
more complex combination between these orbitals could 
exist, but we cannot observe it directly from T(q) with 
6=0. In order to consider other combinations, we set 
9i=9 and change the value of 6. However, even in this 
case, we do not observe any changes in T{q), as observed 
in Fig. 8(b). Namely, the orbital correlation does not 
change due to the rotation in the orbital space, and we 
cannot determine the optimal orbitals. Note that even 
if we optimize 9i at each site, Tiir) is always maximum. 
Thus, we conclude that the states considered in our in- 
vestigations do not have long-range orbital order in the 
PM phase. 

IV. CONCLUSIONS 

In this paper, we have investigated the properties of 
the one-dimensional Hubbard model with three active 
orbitals, with emphasis on electron densities of relevance 
for cobalt oxides. We envision this work as a first step 
toward a numerical accurate study of many-body Hamil- 
tonians for Co-oxides including the Coulombic repulsion. 
Our main result is the identification of two dominant ten- 
dencies in the ground state. For example, at sufficiently 



large Hund's coupling, a tendency toward a fully satu- 
rated ferromagnetic state exists. This state may develop 
long-range charge order in the case when the interorbital 
repulsion U' is large, at density n=5, and for nonzero 
values of the coupling in the direction perpendicular to 
the chains. 

In previous investigations of Co-oxides models, ten- 
dencies toward FM were also discussed (see, for instance, 
Refs ■ fiol and UTI and references therein). Thus, evidence is 
accumulating that magnetic states should be of relevance 
for these materials. This is clearly compatible with ex- 
periments for perovskites cobaltites*^ However, thus far 
only for large Na doping, magnetism has been observed 
experimentally in triangular lattices^ This result may 
arise from the competition with the higher dimensional 
version of the PM state discussed in this work. This state 
has short-range correlations in all channels, and in some 
limits it has an extra SU{4:) symmetry as in two orbital 
models. While this exact symmetry may appear only in 
one dimension and for J=0, remnants may remain under 
more realistic conditions. 

To the extent that our results can be qualitatively ex- 
tended to higher dimensions, the main competition in 
Co-oxide models originates from FM and PM states, with 
long-range and short-range spin and charge order, respec- 
tively. Of course, the effect of geometrical frustration 
could be also an important ingredient to bring about the 
complex spin-charge-orbital structure in the triangular- 
lattice systems. In fact, two of the authors have revealed 
that the spin frustration is suppressed due to the orbital 
ordering in the eg-orbital model on a zigzag chaini^ In 
addition, in the present work we have identified metal- 
insulator transitions with doping away from n=5, while 
the main properties in the spin and charge sectors re- 
main similar as for the integer density n=5. The next 
challenge is to increase the dimensionality of the t2g sys- 
tem toward two dimensions by studying ladders and/or 
zigzag chains. Work is in progress in this direction. 
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